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Applying semi-analytical models of nonideal plasma, we evaluate the behavior 
of the metallic phase in metal-ammonia solutions (MAS). This behavior is mainly 
controlled by the degenerate electron gas, which remains stable down to 5 MPM 
due to high solvent polarizability and strong dielectric screening of solvated ions. 
Comparing the behavior of the metallic state with those of localized solvated elec- 
trons, we have estimated the miscibility gap An for various alkali metals and found 
An(Na)> An(K). It is rather narrow in Rb-NH3 and does not occur in CS-NH3 so- 
lutions, which is in full agreement with the experiments. The case of Li is discussed 
separately. The difference calculated in the excess free energies of the metallic and 
nonmetallic phases is in the order of ksT, yielding a thermally fluctuating mixed 
state at intermediate metal concentrations. It results in a continuous metal-nonmetal 
(MNM) transition above the consolute point T c and a phase separation below T c . 
We propose a criterion for the MNM transition which may be attributed to the line 
of the maximum of compressibility above T c . This line crosses the spinodal one at 
the critical temperature. Finally, we assert that a new electronic phase similar to 
microemulsion should also arise between the spinodal and the binodal lines. 
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I. INTRODUCTION 



Metal-ammonia solutions (MAS) are an example of quantum-classical systems whose 
thermodynamic, dielectric, and optical properties are controlled by electron-electron interac- 
tions. In our recent papers^ we have explained that these interactions are mainly due to the 
dispersion attractions between solvated electrons at metal concentrations in the range of 1-5 
MPM.- These attractions result in two main effects: a phase separation in solutions of light 
alkali metals below a critical temperature, and a dielectric instability of the solution, which 
may be considered as the onset of metallization. Considering the role of the induced dipo- 
lar interactions between solvated electrons, we have restricted ourselves to the nonmetallic 
phase.- We have already given an indication! that the metal-nonmetallic (MNM) transition, 
or insulator-to- metal transition (IMT), are likely driven by a Goldhammer-Herzfeld (GH) 
mechanism,-^ i.e. a polarization catastrophe which has its origin in the dispersion interac- 
tions. We have also indicated that the correlation effects 1 and the solvent polarizability 6 
may be responsible for the unusual properties of the metallic phase and lead to its instability 
at metal concentrations lower than 5 MPM. As explained in Ref.-, this seems to forbid the 
Mott-like scenario for the MNM transition in MAS. That point reactivates the old debates 
of Mott,-&2 Jortner, and Cohen^ 1 ^ 1 ^ on the nature of the MNM transition. Mott assumed 
the phase separation to be a consequence of the Mott-like transition, which is hidden by the 
phase separation below the consolute temperature (T < T c ), and the MNM transition above 
the temperature (T > T c ) should be of the Anderson type due to disorder. On the contrary, 
Jortner and Cohen^ described the system above the consolute point as a microscopically 
inhomogeneous regime in which the concentration fluctuates locally one of the other well- 
defined values M and Mi (M > Mi), and the MNM mechanism does not involve the 
Mott transition. According to them, the concentration Mi « 2 MPM corresponds to the 
nonmetallic blue phase (consisting of solvated electrons as described in Ref .—£■) , whereas the 
concentration Mq k> 9 MPM corresponds to the bronze metallic phase in which electrons are 
delocalized and move freely. Using this hypothesis and a semi-classical theory of percola- 
tion, Jortner and Cohen were able to account for most transport properties observed in the 
intermediate concentration regime above the consolute temperature.— 1 ^ 1 ^ However, in the 
conclusive part of their paper,— they pointed out the four important unsolved questions: 
1) 'what is the origin of this microscopically inhomogeneous state?'; 2) 'what is the origin 
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of the phase separation?'; 3) is there any link between the two phenomena?; and finally 4) 
'what is the nature of the MNM transition?'. We have already partly answered the last 
question in our previous papers^, and we now think that we are able to answer the three 
other ones, at least at a qualitative level at this stage of our theory. This will be the subject 
of the discussion of the present paper, but we are mainly facing a new theoretical situa- 
tion. Our results show that at low temperatures an intermediate concentration range exists 
in which both states, i.e., the nonmetallic one (the blue phase) and the metallic one (the 
bronze phase), are intrinsically unstable.- 1 ^ Such a theoretical possibility has been already 
suggested by one of us^ 3 -^ 4 - 1 ^. As we will discuss in this paper, we believe that this is at the 
origin of the phase separation, and that the existence of an inhomogeneous state above the 
consolute temperature should additionally be clarified. 

There are numerous experimental data on the electronioi&^JLJ&SS and 
structural^i^ 1 ^ 1 ^ 1 ^ properties of the metallic state of MAS, but a theoretical treat- 
ment of concentrated MAS is restricted to a few papers only— iiii2iii2L2&i2ii The difficulty of 
such a treatment is twofold: first, the electron density in the metallic phase is still rather 
low so that strong electronic correlation effects occur, and, secondly, the molecular nature 
of the solvent as well as a the sufficient role of the solvent polarizability complicate the 
theoretical study of these correlations in the metallic phase. The aim of the present work 
is to develop a statistical model treating the excess electrons in MAS under conditions 
corresponding to the metallic phase. This phase represents a three-component mixture 
consisting of delocalized electrons, cations, and polarizable solvent molecules. In general, a 
detailed information about electronic and thermodynamic properties of the system may be 
obtained by quantum molecular simulations based on the Car-Parrinello method.— i22a2L22 
However, the computational costs and sensitivity of such calculations to the choice of 
the interaction potentials (an explicit account of cations, polarizability of solvent and so 
on) restrict the application of the method and cannot be used efficiently for a complete 
understanding of the phase behavior of MAS. An alternative way may be an employment 
of the integral equations theory based on quantum hypernetted chain closure 33 , but 
the current status of such an approach is limited only to simple metals 34 or electrons 
localized in polar liquids 3 ^ 1 ^ and ionic liquids.— 1 ^ Indeed, we do not know any applications 
extended to molecular metallic fluids. In order to avoid the complexity related to numerical 
implementations, we consider in the present paper that the metallic phase is well described 
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by an effective two-component plasma (TCP), where the role of the solvent is restricted 
to its screening effects. We take into account different frequency scales of this screening 
as we have done for the nonmetallic phase,-" 2 - i.e. that the screening of electron-electron 
interactions is due to the electronic polarizability of solvent, whereas ion- ion interactions 
are additionally screened by the permanent dipoles of the solvent molecules. As a result, 
we use different dielectric constants to treat the Coulomb interactions properly, i.e. a 
high-frequency dielectric constant for the electron-electron interactions and a low-frequency 
dielectric constant for ion-ion interactions respectively. By using available analytical 
expressions for the free-energy of the TCP— we will derive the main thermodynamic and 
dielectric characteristics of the metallic phase, its pressure and chemical potential, and we 
will evaluate the locus of the critical lines (spinodal and coexistence line) which correspond 
to the instability of the metallic phase. Then, by comparing the thermodynamic behavior 
of the metallic phase with that of the nonmetallic phase,- 12 we will give evidence that an 
inhomogeneous electronic microstructure must arise at intermediate metal concentrations in 
MAS. Although the detailed study of the inhomogeneous microstructure will be the subject 
of our next paper,— we will assert in this work that the thermodynamic and dielectric 
peculiarities of the inhomogeneous electronic state govern the behavior of MAS in this 
concentration range. Eventually, by treating the inhomogeneous microstructure within the 
methods of simple liquids, we will propose a macroscopic criterion for the MNM transition, 
and sketch the complete phase diagram of MAS. Atomic units are used throughout. 



II. MODEL 



A. General outline. 



The high-density phases of MAS are composed of delocalized interacting electrons scat- 
tered by ammonia molecules and solvated ions. As explained above, the solvent molecules 
are considered to result in the screening of the ions and the electrons only. The ion-ion 
interactions will be screened by the low frequency dielectric constant e/, while the electron- 
electron interactions will be screened only by the high frequency dielectric constant e^. Next 
we treat the system as a TCP, the first component being the interacting electron gas in a 
jellium, and the second component being a classical plasma of ions in a jellium of electrons. 
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Finally, we add an interaction term which represents a correction to the Madelung energy 
of the point ions in the jellium of electrons, i.e. a pseudo-potential correction due to the 
short-range scattering of the free electrons by valence core electrons of the ions. Notice that 
these ions are themselves solvated by ammonia molecules. This defines two different effec- 
tive radii for the ions. The first one is an effective classical radius of the solvated ions, that 
we call r V( iW) which will be used to account for the classical short-range ion- ion interactions 
through the use of a packing factor. The second effective radius concerns the electron-ion 
interactions and it is related to the core ion radius, and we call it ?v It includes the effect of 
the solvation shell of a particular ion. For alkali metals such as Na, K, Rb, Cs, the effective 
core radii are only slightly different from those determined by the pseudo-potential method 42 
for non solvated ions. However, the case of Li + is special, since it is well established 2 ^ that 
owing to their small size, the Li + ions occupy tetrahedral vacancies formed exactly by four 
ammonia molecules. From this point of view, Li + ions are strongly bound to form complexes 
Li(NH3)4~, which are dominant ionic species in the solution. Consequently the core radius 
of solvated Li + is substantially larger than in the case of liquid Li metal. 

Thus, our model contains four parameters ei, eh, r v aw, r ion, in addition to the Wigner-Seitz 
parameter r s = fr 2 /me 2 (47m/3) 1//3 , related to the metal density (or the density of electrons, 
since the alkali metals are monovalent). Finally, the excess free energy per electron of the 
metallic phase is written as the sum of three terms: 

f m (n, T) = e DEG (n) + focp(n, T) + f e i(n), (1) 

where eoEG{ n ) is the f ree energy per electron of an interacting electron gas with density n, 
and focp{n,T) is that of a classical one-component plasma (OCP) of ions in a jellium of 
electrons. Finally, the last term f e i(n) represents the electron- ion interactions and will be 
treated in the framework of the pseudo-potential theory 4 ^. Notice that Ashcroft proposed a 
similar approach to the metallic phase,— except that he didn't used the different dielectric 
constants in his model. 

B. Free-energy of the interacting electron gas 

The electron gas may be parameterized by two dimensionless parameters 3 ^: the elec- 
tron coupling constant T e = / 9e 2 (47rn/3) 1 ' 3 /e/ l , and the reduced temperature 6 = 
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2m/h 2 f3(3ir 2 n) 2 / 3 , where (3 is the inverse temperature, m and e are the electron mass and 
charge, respectively. However one can easily check that T e ~ 100 ^> 1 and ~ 0.05 <C 1 
under conditions corresponding to the metallic phase in MAS, and therefore the electron 
gas is degenerate and we can treat it at T = OK. Its free energy may be parameterized by 
the single parameter, i.e. the effective Wigner-Seitz radius defined as r* = r s /t h . Applying 
the conventional methods for the electron gas, we have: 

(t 2 h Ti 2 \ . . 1.105 0.458 

The first two terms are obtained by the usual Hartee-Fock approximation, and the last term 
is the correlation energy (beyond the Hartree-Fock approximation). In order to treat this 
correlation energy e cor (r*), which is a smooth function of r*, we have employed the result of 
the local density approximation parameterized in Ref.— by fitting quantum simulations: 

~ 70 (3) 



^cory s 1 -I In ; 

1 + 7ir* i ' 2 + 7 2 r* 

where 70 = —0.1423, 71 = 1.0529, and 72 = 0.3334 are numerical parameters obtained from 
ref.— 3 . As we noticed in Ref.-, the degenerate electron gas gives the main contribution to 
the total excess free energy f m (n,T) of the metallic phase. 



C. The classical one-component plasma 

The mean distance between ions exceeds 9 A at metal concentrations of about 10 MPM, 
so short-range details of the ion-ion interactions can be ignored, and we may treat the ions 
as charged hard-spheres in the uniform jellium background. Thus, the excess free energy 
of ions is controlled by the two dimensionless parameters: the dimensionless ion coupling 
constant Tj = / 9e 2 (47rn/3) 1/ ' 3 /e;, and the packing factor rj = ima 3 /6 in which a = 2r v aw is 
the classical hard-sphere ion diameter. Simple evaluations show that the first parameter is 
of about 10 in the metallic Na-NH 3 , and that the second one does not exceed 0.05 under 
the same conditions. Therefore the free energy focp per ion may be written as 

Pfocp = MnA 3 ) - 1] + f hs (v) + fci^i), (4) 

where the first term is the ideal contribution, A$ = (2nf3 / 'Mi) 1 / 2 is the de Broglie length and 
Mi is the ion mass. The second term in (T4J) is the hard-sphere contribution, and the last 
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term is the contribution due to Coulomb interaction between ions. Analytical expressions 
for these contributions are well-known. The hard-sphere part is expressed as^: 

Pfhsiv) = (1 _^ a ■ (5) 

According to Ref.— , in which simulation data were fitted for various Tj in the range 1 < 
Tj < 160, the electrostatic contribution can be approximately written as: 

Pfc(Ti) = -0.897^ + 3.620r- /4 - 0.758iy 1/4 - 0.815 In T { - 2.58. (6) 

The first term in this expression represents the Madelung energy and gives the most im- 
portant contribution. The remaining terms are temperature-dependent corrections due to 
thermal motion of the ions. 



D. Electron-ion interaction 

The main difficulty is to evaluate f e i(n), because the bare ion-electron interaction is 
modified in the solution by the polarizability of the solvent. Moreover, the Madelung term 
in eq. concerns the point ions. The correction to this approximation is well treated 
by the pseudo-potential model for simple metals^. In that point we follow the Ashcroft 
approach^ and adapt the pseudo-potential model^ to our case. As a result, we write the 
electron-ion contribution (expressed in effective atomic units me^/effi 2 ) as: 

(§) /-<»> - §?• e) 

where r* = r^/e^ and the ion-core parameter r, is related to the atomic number of the ion 
as discussed above. Notice that in our previous paper,- we took a slightly different notation 
and wrote f e i(n) = din, so that aj = Qnr 2 to make the exact link with the present paper. 

III. THERMODYNAMIC AL PROPERTIES 

Once we know the expression of the excess free energy f m per electron, we can deduce 
the excess pressure Ap, the excess chemical potential A/x, and the excess compressibility k 
as determined by the usual formulas: 

AP = H ^ ^ = fm + n—, k =n—. (8) 
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As discussed above, it may be checked that the degenerate electron gas (DEG) gives the 
main contribution to f m (n,T) at low temperatures, whereas the ionic contribution is only 
a correction. Our TCP model of the metallic phase of MAS depends on temperature and 
metal concentration and can be characterized by the four parameters: r v dw, r i, Q> and e^. 
As we will see, the value of the ionic mass M, in (j4j) does not change the values of the density 
at which the metallic state becomes unstable. 

Earlier on2, we have determined the low-density spinodal line n s (T) above which the 
solvated electrons become thermodynamically unstable. Now we evaluate the high-density 
counterpart of the spinodal line n C 2(T) corresponding to the zero derivative of the excess 
chemical potentials or the excess pressure: 



Since the main part of the excess free energy in the metallic phase is coming from the 
degenerate electron gas, the high-density spinodal n C 2(T) will depend weakly on temperature. 

For any second-order phase transition there are two curves: the spinodal and the bin- 
odal ones, the last curve corresponds to the liquid-liquid coexistence. As we explained in 
Ref.-, the low- and high-density spinodals are given by n s (T) and n C 2(T), respectively. In 
principle, to calculate the coexistence curve, we need a complete knowledge of the excess 
free energy f(n,T) in the whole range of density, which is beyond the current status of our 
theory. However, we may use the following argument to estimate the high-density part of 
the coexistence curve: the excess pressures of the low-density and high-density phases must 
be equal (this is the coexistence condition). But we know that the excess pressure of the low 
density phase is roughly close to zero (more exactly it is in the order of nksT), 2 whereas 
the excess pressure of the metallic phase strongly varies with n. Therefore, we assert that 
the coexistence curve n p (T) at high density will be roughly given by the condition: 



The critical concentration n C 2(T) defines the low-density limit of the absolute stability 
of the metallic phase, whereas the critical concentration n p (T) is close to the line of liquid- 
liquid coexistence. The metallic phase cannot exist below n C 2(T), and in the range n C 2(T) < 
n < n p (T) the solution is not homogeneous (the pressure becomes negative) and consists 
of domains of the metallic phase together with nonmetallic islands being nuclei of the low 



dAji/dn = dAp/dn = 0. 



(9) 



Ap(n p ) = 0. 



(10) 
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density phase, i.e. the solvated electron phase. The uniform metallic phase remains stable 
for densities n > n p (T). 

IV. THE USE OF EXPERIMENTAL DATA TO FIT THE MODEL 

PARAMETERS 

Our input phenomeno logical parameters are r i; r v dw, £h, an d ej. For the pure ammonia 
solvent, we take = 1.756, e s (T = -70°C) = 25, and de s /dT = -0.1 K" 1 .^ We use 
the data on the ammonia density,— namely, riNH :i — 0.0255A~ 3 at T = — 70°C. In order to 
take into account phenomeno logically the influence of free electrons on the high-frequency 
dielectric constant, we apply: 

eh(n) = €oo - Ac, (11) 

where c = n/{n + n^H-i) is the relative fraction of the metal in the solution, and A is 
a phenomenological parameter. In general, the low-frequency dielectric constant e\ may 
deviate also from the bulk value e s due to saturation of orientational polarization in the 
vicinity of ions. Simple estimates show that the low-frequency constant can decrease down 
to 9. 48 Therefore, different values for ei could also be tested. The core-radii and van-der- 
Waals radii can be extracted from the literature, which we do in our first parameterization 
(model 1). However, as we discussed above, this does not take properly into account the 
fact that the ions are solvated. Then we propose a second procedure to evaluate these radii, 
that we call model 2. 

Model 1. The ion-core radii rj are obtained from the data on simple metals.— Although 
we do not find suitable estimates of van-der-Waals radii r V( j,w for ions solvated in ammonia, 
we have used those of hydrated ions, assuming that the solvation properties of water are not 
so different from those of ammonia. We take the values for r v dw from Ref.—. That gives 
the values reported in Table 1 (at the lines: model 1). Finally, in this first model, we apply 
the phenomenological parameter A = 1 — = 0.756, that takes into account the volume 
fraction occupied by the metal as in Ref.—. 

Model 2. For this model, the phenomenological constant A is used as a fitting parameter 
of experimental data on plasmon excitation measured in Li-NH3.— It is experimentally 
found €h = 1.44 at 20 MPM so that A = 1.58. As discussed above, the case of Li raises a 
question, since the core radii obtained from Ref.— do not take into account the solvation shell 
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of the ion. We use as a parameter to fit the experimental data on the phase coexistence, 
which is known for Na-NH 3 , Li-NH 3 , and K-NH 3 solutions.— We fit these data by (ITU1) . 
and obtain the values for the effective ion-core radius. A similar method is often used to 
fit ion-core radii for simple liquid metals.— Comparing T{ to r v dw for model 1, we find a 
proportionality constant of about 3/2 for Na-NH 3 . Then we use the same proportionality 
constant to deduce r v dw from for Li, K, and Rb, which gives the values reported in Table 
1 (at lines: model 2). 

Table 1. Parameters of the metal ions and the calculated critical densities for models 1 and 2. 





model 1 


model 2 


ion 


r,(A)i2 


w(A)i2 


n s (MPM) 


n c2 (MPM) 






n s (MPM) 


n c2 (MPM) 


Li 


0.761 


1.33 


1.1 


7.1 


1.25 


1.88 


1.72 


3.8 


Na 


0.939 


1.68 


1.32 


4.92 


1.12 


1.68 


1.32 


5.2 


K 


1.234 


2.02 


2.04 


3.0 


1.31 


1.97 


1.96 


3.4 


Rb 


1.323 


2.16 


2.76 


2.65 


1.33 


2.0 


1.96 


3.2 


Cs 


1.449 


2.54 


5.0 


2.15 











V. ANALYSIS 



A. Pure metallic phase 

Using the formulas described above, we have calculated the excess chemical potential 
and the excess pressure for the metallic phase of Na-NH 3 solutions. Figure 1 shows the con- 
centration dependencies of these thermodynamic characteristics at various temperatures for 
model 1. As can be seen, the chemical potential has the minimum at n c2 ~ 5 MPM and the 
pressure is negative below n cp = 9 MPM at low temperatures. Therefore, the metallic phase 
is not stable at densities lower than n c2 . This instability results from Coulomb correlations 
between delocalized electrons, because the main contribution to this compressibility is due 
to the DEC It is well known (see f.e.— ) that the compressibility of a pure DEG in jellium 
(without polarizable medium), becomes negative at r s > 5.24 and its pressure is negative for 
r s > 4.18, which would correspond to n c2 = 30 MPM and n p = 46 MPM, respectively. These 
estimates are far from the experimental values in MAS, which clearly indicates that the high 
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solvent polarizability and strong dielectric screening of ions stabilize the DEG in MAS to 
lower metal concentrations. The temperature influence on the stability is due to the ionic 
contribution, and it is as small as it should be for a quantum phase transition. The critical 
concentration n C 2 varies only from 5.5 to 4.5 MPM as the temperature changes by 100° C. 
The critical concentrations obtained are sensitive to the choice of the input parameters rj, 
e s , and A. Figures 2-4 show the variations of these critical concentrations with respect to 
these parameters. A decreasing contribution of the solvent polarizability tends mainly to 
destabilize the metallic phase, since n c2 oc e^ 3 (n), and an increasing ion-core radius promotes 
the phase stability at lower densities. 

In our theory, we can estimate the miscibility gap by considering the range between the 
critical lines obtained at low density n s (T) (our previous paper-) and the present calculations 
of n c2 (T). Both the lines n s (T) and n C 2(T) give the locus of the low- and high-density 
spinodals, while the range n s < n < n c2 correspond to the miscibility gap. The results 
calculated at T = — 70° C are listed in Table 1. We find that the gap decreases monotonically 
as the ion size rises as it happens for simple metals and the gap disappears for heavy ions 
like Cs, which is experimentally confirmed.^ Although these qualitative predictions provide 
a correct trend for the miscibility gap versus the ion size, the estimates based on the ion-core 
parameter derived from the theory of simple metals^ (our model 1) do not yield accurate 
evaluations of the phase separation range in case of Li-NH 3 . The reason of this discrepancy 
is a peculiarity of local microstructure around the solvated ions, as we discussed above. 
Namely, the size of the Li ion is small, and the coordination number for the solvated Li ions 
is small too in comparison to other ions. As a result, the delocalized electrons scatter on the 
solvated Li-NH 3 complexes rather than on the cores of Li ions, and the effective ion-core size 
deviates sufficiently from the value obtained for simple metals. This drawback disappears 
when we apply model 2 and use the data on rj fitting the zero-pressure condition (flUj) at the 
coexistence line. The critical concentrations n s and n C 2 calculated in this manner correspond 
much better to the experimental data (see Table 1); in particular, they yield the correct trend 
in the range An = n C 2 — n s of the phase separation, namely, An(Na)> An(Li)> An(K). 
Both models 1 and 2 also indicate the absence of the phase separation in Cs-NH 3 solutions, 
while our results for Rb-NH 3 solutions exhibit a possibility a phase separation for model 2 
and absence for model 1. The experimental evidence of such separation is not clear. There 
is no visible phase separation like in the Li-NH 3 solutions,— although the measurements 
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of conductivity^ and X-ray scattering^ indicate large fluctuations at n c rj 4 MPM and 
at temperatures close to T c = 197°K. Those fluctuations can be interpreted as an indirect 
evidence of a phase separation for RD-NH3 solutions in this range.— ^ 

We have also evaluated the compressibility k as a function of the metal concentration n 
in the framework of model 2. The inverse compressibility k~~ x calculated is shown in Fig. 5 
together with the experimental data derived from the measurements of plasmon excitations^ 
obtained for the metallic Na-NH 3 and Li-NH 3 solutions above the consolute point T c . As 
can be seen, there is a qualitative agreement between the calculated and the experimental 
data. Both indicate that the inverse compressibility decreases as the metal concentration 
decreases. However, the theoretical evaluation gives a much higher compressibility than the 
experimental measurements at intermediate densities. The deviations are more pronounced 
as the concentration decreases. In our opinion, it is an indirect evidence of the formation 
of an inhomogeneous electronic state consisting of a microscopic mixture of delocalized and 
localized electrons, as will be discussed in Sec. IV Bl 

Finally, with the use of the calculated compressibility we can evaluate the low-frequency 
dielectric response of the metallic state at small wave- vectors (k —> 0), since we have the 
limiting relation for this function, as it takes place for the usual electron gas: 

e ( w ^0,*^0)=^n) + ^*p, (12) 

where qxF is the Thomas-Fermi screening wave vector equal to 2(3n/n) 1 ^ 6 /a^ (ao being the 
Bohr radius), and K/ ree = (3 / ir^ri) 1 ' 3 / n is the compressibility of the ideal Fermi gas. The 
function s~ l {k) = l/e(u = 0, k — > 0) is depicted in Fig. 6 at various metal concentrations. 
This function becomes negative at concentrations below n C 2(T), which confirms that the 
region of phase separation is indeed a region where the static dielectric constant must be 
negative.— 1 ^ 

The experimental phase diagram of sodium- ammonia solutions^ 1 ^ 1 ^ together with the 
different calculated critical lines corresponding to various instabilities of the nonmetallic 
and the metallic phases are depicted in Fig. 7. The low-density and the high-density spin- 
odals correspond to the critical lines n s (T) and n C 2(T), respectively, and the line n c \{T) 
of polarization catastrophe (see paper-) gives the onset of the MNM transition. Although 
our calculations of critical lines roughly correspond to the experimental situation, the ex- 
perimental coexistence line coincides with the theoretical zero-pressure line of the metallic 
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phase only at low temperatures and deviates significantly from it as the temperature rises. 
The situation is similar for the calculated low-density and high-density parts of the spinodal 
curve, they do not cross as they should at the consolute point T = T c . Therefore, at this 
stage of the theory, our model of the uniform metallic state and the homogeneous phase of 
solvated electrons phase cannot explain all the peculiarities of the phase behavior of MAS, 
although it indicates the main features of this behavior, namely, the existence of a miscibility 
gap giving rise to a phase separation. 



B. Thermally fluctuating inhomogeneous state 

To qualitatively understand what happens in the intermediate region, we have calculated 
the difference f nm — f m between the excess free energies of the nonmetallic and the metallic 
phases. Following our previous study,- the free energy per electron in the nonmetallic state 
can be expressed as: 

/ -™ (n) = h " t " t)^ + 4lCsAr - + — sr - ^ + (13) 
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e 

„3 



ln(nAf 2 Af ) - 1 + ^ ~ ^ 

\ -L ' tri m , 



nm ) 
2 



2 Y C a nr* 



e s[°" + 7 x ] 3irn(3 e 



2 

DO 



Vnm ) 

where Cy, Cs, and C a are numerical coefficients that we have already discussed in Ref.- 
and related to the temperature- and concentration- dependencies of the radius r e (n,T) of 
the solvated electrons, A e = (2nP/M e ) 1 ^ 2 is the de Broglie length for localized electrons, M e 
is their effective (classical) mass, A is the surface tension, and r) nm = im(r e + r v ow) 3 /3 is an 
effective packing factor, that we take equal to the mean value between the solvated electron 
and ion diameters. Eventually, 7 = ([1 + (327m/e s ) 1 / 2 (r e + r v r>w)} 1 ^ 2 — l)/2(r e + r v ow) is 
the inverse screening length. The first row of the expression yields the solvation free energy 
of noninteracting electrons, and the second row results from electron-electron interactions 
and includes short-range, electrostatic, and dispersion contributions, respectively. All these 
contributions are only corrections to the free-energy of the electron solvation energy. 2 

To calculate f nm , we have used our model described in Ref.-, namely: = 1.5, C r ~ 1.25, 
C s = 1 and C v = 1.75, A = 40 dyn/cm, and we apply C a = 0.14 and A e = 0.27A. The 
difference Sf(n) = f n m{n) — / m (n) is shown in Fig. 8a. The main point to underline is that 
their difference remains in the order of only a few fcgT along the range of concentration 1-10 
MPM, although the values of the free energies of the metallic and non-metallic state are 
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both in the order of -0.7 eV. We should emphasize that this is not by chance. By carefully 
examining (and simplifying) the energy of both the metallic and non-metallic states, we can 
indicate that the main part of the energies coming from the electronic part is roughly given 
by: 

f ~ 112 ( 1 I^M f ~ H °- 46 i e -or{r s /e^) 

Jnm ~ 9 I I ; Jm ~ 9 + 9 ; I "J 

where / m stands for the metallic phase, and f nm for the nonmetallic one. That clearly 
indicates a crossing between the two energies when r s tends to r e , which is in quantitative 
agreement with the results shown in Fig. 8a. 

Consequently, the thermal fluctuations between these electronic states must play a cen- 
tral role at intermediate densities, in particular to allow the closure of the miscibility gap at 
the critical temperature T c by a mixing entropy effect. The detailed calculations of the com- 
plete phase diagram, including the calculation of T c will be reported in our future paper.— 
However, at a qualitative level, we may say that above the critical temperature, the system 
should be a microscopic mixture of both states, which is highly thermally fluctuating and 
roughly described by the relative fractions n m = (1 + exp[/3(/ m (n) — fnmij 1 ))]) -1 an d 1 — n m 
of electrons in the metallic and nonmetallic state respectively, as a function of the metal 
concentration. The calculated values of n m are also shown in Fig. 8b for Na-NH 3 solutions 
at T = 240° K. They indicate that our initial criterion of metallization, i.e. the polarization 
catastrophe above n c i(T)^ is spread out above T c and becomes a progressive MNM transi- 
tion, which should take place around the density n c i(T). Another peculiarity of the MNM 
transition is the compressibility of the mixed state consisting of localized solvated electrons 
as well as delocalized free electrons. The line of maximum of compressibility (or minimum 
of the inverse compressibility) may well characterize the MNM transition. This is sketched 
on Fig. 9 by the total inverse compressibility of the mixture. The total compressibility is 
well given by the sum of the inverse compressibility for each species (solvated electrons and 
delocalized ones) at the respective densities 1 — n m and n m . At the same time, the dielectric 
constant remains finite at the line of the minimum inverse compressibility, although it rises 
sharply in the nonmetallic phase close to this line, as discussed in our previous paper-. It 
is worth noting that the mixed state qualitatively described here, is of a new kind in con- 
densed matter (to our knowledge), and really results from the competition, driven by strong 
thermal fluctuations, between the delocalized and the self-trapped quantum states. 
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The situation is more simple below T c . There is a finite range of densities between the 
spinodal and binodal curves at each side of the phase diagram, where the minority phase may 
nucleate in the majority phase to form some kind of stable electronic microemulsion. This 
microemulsion phase would be characterized by a large variety of aggregates as reported 
in numerical simulations 2 - 9 . At the same time, the solution is phase separated below the 
spinodal line to prevent negative compressibility. The Fig. 10 summarizes all these effects. 
The MNM transition may be assigned to the line of the minimum inverse compressibility 
which crosses the spinodal and coexistence lines at the critical point n c (T c ). A more detailed 
theory of all these effects will be proposed in our future publication as well as more precise 
calculations of the spinodal and coexistence lines^. We have only qualitatively discussed in 
this paragraph the main features of the phase diagram resulting from our model. 



VI. CONCLUSION 



Thus, using the methods of non ideal plasma, we have evaluated the behavior of the pure 
metallic phase in MAS. We have asserted that this behavior is mainly controlled by the 
state of degenerate electron gas. Due to high polarizability of ammonia and high dielectric 
screening of ions solvated in MAS, the gas remains stable up to values of r s of about 11, 
which is quite different from simple metals where available values of r s do not exceed 6. 
The pure metallic phase is unstable at metal concentrations lower than 5 MPM. Comparing 
this critical concentration n c2 with the critical concentration n s (T) corresponding to the 
van-der-Waals instability of solvated electrons, we have evaluated the range An of the phase 
separation and found it to be governed by the ion-core size, the range decreasing as the ion- 
core size rises. By evaluating the ion-core size with the use of the zero-pressure condition 
applied to the coexistence line, we have obtained the correct trend in the range of the phase 
separation, namely, An(Na)> An(Li)> An(K). The phase separation does not occur in 
Cs-NH 3 solutions and is rather narrow in the case of Rb-NH 3 . 

Comparing the calculated free energy of pure metallic and nonmetallic phases, we have 
asserted the difference between the energies to be small due to minor differences in the en- 
ergy of polaron formation and that of electron gas at the stability boundary. This leads to a 
strong influence of thermal fluctuations on the nature of the MNM transition, which becomes 
continuous above the critical temperature T c . A thermally fluctuating inhomogeneous elec- 
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tronic mixture arises. It modifies substantially the criterion for the MNM transition, because 
this transition should be attributed to the line of the minimum inverse compressibility or 
to the maximum scattering factor. This line is quite close to our previous estimates n c x(T) 
corresponding to the polarization catastrophe and divergent dielectric constant, 2 although 
the dielectric constant remains finite at temperatures exceeding T c . 

At the same time, the dielectric constant and the compressibility must be negative below 
the spinodal line. To avoid such negativity, the system macroscopically separates into the 
metallic and nonmetallic phases below the spinodal line. A new electronic state arises in 
the range restricted by the spinodal and the binodal lines. This state can be characterized 
as an electronic microemulsion. The locus of the critical point n c and the spinodal line is 
determined by the thermodynamic behavior of microemulsion and will be the subject of 
our next study— Concerning the criterion for the MNM transition, we should note that it 
is quite different from the one for usual semiconductors^ and indicates a sufficient role of 
thermal fluctuations in the mechanism of MNM transitions. On the other hand, the line 
of the minimum inverse compressibility always exists for the second order transitions in 
simple^ and molecular^ liquids as well as for self-assembled networks.— This line crosses 
the spinodal line at the critical point. That is the reason why the MNM transition couples 
with the phase separation in MAS. 
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FIG. 1: The excess pressure Ap(n,T)/n (a) and the excess chemical potential Afi(n) (b) versus 
metal concentration n at various temperatures for the metallic phase of Na-NH3 solution (model 
1). 




FIG. 2: Critical concentrations n C 2 and n p versus the parameter A for the metallic phase of Na-NH3 
at T = -70°C (model 1). 
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FIG. 3: Critical concentrations n C 2 and n p versus low- frequency dielectric constant e s for the 
metallic phase of Na-NH 3 at T = -70°C (model 1). 




FIG. 4: Influence of the ion-core radius on the critical concentrations n C 2 and n p in the case of the 
metallic phase of Na-NH 3 at T = -70°C (model 1). 
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FIG. 5: The concentration dependence of the inverse compressibility k _1 in MAS at T = 240° 
K. The solid curve indicates our calculations in the case of L1-NH3 (model 2), and the diamonds 
are the experimental data from2£, the dashed curve and the squares are the same in the case of 
Na-NH3 solutions. 
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FIG. 6: The inverse static dielectric function e~ 1 (k,co = 0) at various metal concentrations in 
L1-NH3 solutions at T = 240° K (model 2). The values of metal concentration n are indicated at 
the corresponding lines. 
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FIG. 7: The phase diagram of the Na-NH3 solution and the lines of instabilities calculated by 
model 2. The experimental dat a 51 ' 55 ' 56 are indicated by symbols, and their cubic interpolation 
is indicated by the thick curve. The dotted line corresponds to the polarization catastrophe n c \ 
considered as the onset of metallization. The thin solid line is the low-density spinodal n s , the 
dashed one indicates the high-density spinodal n C 2, and the dashed-dotted curve yields the zero- 
pressure line of the pure metallic phase. The bottom of the figure corresponds to the solidification 
of ammonia. 
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FIG. 8: The dimensionless difference /3(f m — f nm ) in the free energies of metallic and nonmetallic 
phases (a) and the relative fraction of electrons in delocalized states (b) versus the metal concen- 
tration in the Na-NH3 solution at T = 240° K, all other parameters of the calculations correspond 
to Fie. 1. 
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FIG. 9: Schematic sketch of the influence of thermal fluctuations on the concentration dependencies 
of the inverse compressibility above (a) and below (b) critical temperature. The dashed-dotted 
line corresponds to localized electrons, the dashed line represents the inverse compressibility for 
delocalized electrons. The solid curves correspond to the total inverse compressibility obtained as 
a sum of the contributions of delocalized and localized electrons, proportionally to their respective 
densities. 
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FIG. 10: Schematic sketch of the phase behavior of Na-NH3 solutions. The dashed line corresponds 
to the MNM transition, which is attributed to the line of the minimum inverse compressibility. The 
solid curve is the spinodal line. The dashed-dotted curve corresponds to the line of liquid-liquid 
coexistence. 



